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We demonstrate by numerical flux calculations that neutrino beams producing the observed 
highest energy cosmic rays by weak interactions with the relic neutrino background require a non- 
uniform distribution of sources. Such sources have to accelerate protons at least up to 10 23 eV, 
have to be opaque to their primary protons, and should emit the secondary photons unavoidably 
produced together with the neutrinos only in the sub-MeV region to avoid conflict with the diffuse 
7 —ray background measured by the EGRET experiment. Even if such a source class exists, the 
resulting large uncertainties in the parameters involved in this scenario does currently not allow to 
extract any meaningful information on absolute neutrino masses. 


I. INTRODUCTION 


In acceleration scenarios ultra high energy cosmic rays 
(UHECRs) with energies above 10 18 eV are assumed to 
be protons accelerated in powerful astrophysical sources. 
During their propagation, for energies above > 50 EeV 
(1 EeV = 10 18 eV) they lose energy by pion production 
and pair production (protons only) on the cosmic mi¬ 
crowave background (CMB). For sources further away 
than a few dozen Mpc this would predict a break in 
the cosmic ray flux known as Greisen-Zatsepin-Kuzmin 
(GZK) cutoff [ji] , around 50 EeV. This break has not been 
observed by experiments such as Fly’s Eye Q, Haverah 
Park 0 , Yakutsk Q and AGASA 0, which instead show 
an extension beyond the expected GZK cutoff and events 
above 100 EeV. However the new experiment HiRes 10 ] 
currently seems to see a cutoff in the monocular data [|7j. 
Taking into account that all old experiments except per¬ 
haps AGASA do not have sufficient statistics in the high¬ 
est energy region to settle the question, the existence of 
a possible cutoff remains unclear at the moment. The 
apparent absence of a cutoff especially in the AGASA 
data has in recent years triggered many theoretical ex¬ 
planations ranging from conventional acceleration in as¬ 
trophysical sources to models invoking new physics such 
as the top-down scenarios in which energetic particles are 
produced in the decay of massive relics from the early 
Universe This enigma has also fostered the develop¬ 
ment of large new detectors of ultra-high energy cosmic 
rays which will increase very significantly the statistics 
at the highest energies |J. 

In bottom-up scenarios of UHECR origin, in which 
protons are accelerated in powerful astrophysical objects 
such as hot spots of radio galaxies and active galactic 
nuclei |RJ , one would expect to see the source in the di¬ 
rection of arrival of UHECRs, but above the GZK cutoff 


in general no suitable candidates have been found within 
the typical energy loss distance of a few tens of Mpc 
for the known electromagnetically or strongly interacting 
particles . Even assuming significant deflection by 

large scale extragalactic magnetic fields requires at least 
several sources 13] whose locations have not been iden¬ 
tified yet. 

Moreover, recent observations of small scale cluster¬ 
ing by the AGASA experiment L 0 suggest that sources 
of UHECR are point-like |]l^,^6[ . This fact together with 
the lack of nearby sources favors the possibility of sources 
much further away than 100 Mpc, at redshifts of or¬ 
der unity. An additional motivation for this possibility 
comes from recently reported possible correlations of the 
arrival directions of observed UHECR above ~ 50 EeV 
with certain classes of sources such as compact radio 
galaxies |n]] or BL Lacertae objects [[Is]. In the latter 
case it is still possible that the sources are located at 
moderate distances z ~ 0.1. In this case photons with 
extremely high energies E > 10 23 eV can propagate sev¬ 
eral hundred Mpc (constantly loosing energy) and can 
create secondary photons inside the GZK volume [ |l9f . 
However, this model requires both extreme energies of 
primary photons and extremely small extra galactic mag¬ 
netic fields (EGMFs) B < 10~ 12 G. Moreover, if a corre¬ 
lation with any source at redshift z > 0.2 is found, this 
model will be ruled out. 

If sources of the highest energy cosmic rays are in¬ 
deed at cosmological distances £ ~ 1, the only known 
mechanism not involving new physics except for neu¬ 
trino masses assumes neutrinos as messenger particles: 
Charged particles accelerated in such sources give rise to 
a secondary neutrino beam which can propagate essen¬ 
tially unattenuated. If this neutrino beam is sufficiently 
strong it can produce the observed UHECRs within 100 
Mpc by electroweak (EW) interactions with the relic neu¬ 
trino background |2(J. Specifically, if the relic neutrinos 
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have a mass m„, Z-bosons, whose decay products can 
contribute to the UHECR flux, can be resonantly pro¬ 
duced by ultra high energy (UHE) neutrinos of energy 
E v ~ M 2 z /{2m v ) ~ 4.2 x 10 21 eV (eV/m„). 

However, this “Z-burst” mechanism is severely con¬ 
strained by at least two types of observational data: 
First, there are upper limits on the UHE neutrino flux, 
based on the non-observation of horizontal air showers 
by the old Fly’s Eye experiment [^TJ or by the AGASA 
experiment j 22 [ and from the non-observation of radio 
pulses that would be emitted from the showers initiated 
by the UHE neutrinos on the moons rim [^3|. Second, 
even if the sources exclusively emit neutrinos, the EW in¬ 
teractions also produce photons and electrons which ini¬ 
tiate an electromagnetic (EM) cascade which transfers 
the injected energy down to below the pair production 
threshold for photons on the CMB || . The cascade thus 
gives rise to a diffuse photon flux in the GeV range which 
is constrained by the flux observed by the EGRET in¬ 
strument on board the Compton 7 —ray observatory [ p4| . 
Reproducing the observed UHECR flux by the Z-burst 
mechanism under these two constraints has been shown 
to in general require local relic neutrino over-densities 
in order to increase the local UHECR flux. These over¬ 
densities turn out to be much higher than values 2-3 
which would be expected from the over-density in the 
local supercluster fl25|| . 

In order to avoid this difficulty one can suppose that 
the Z-burst mechanism is responsible only for part of the 
UHECR flux |2(|. In this case, one can reduce both pri¬ 
mary neutrino and secondary photon fluxes and obey all 
existing limits. However, the price for this is to explain 
only a part of the UHECR events by the Z-burst mech¬ 
anism and the necessity for a second source mechanism 
for UHECRs. 

Furthermore, Ref. ^ 6 | claims that already the present 
data provides possible evidence for the relic neutrino 
background and starts to constrain the absolute neutrino 
mass, a possibility that has recently been discussed in 
principle in Ref. p9| . This claim is based on tuning many 
unknown parameters such as the value of the EGMF, the 
universal radio background (URB) which governs pair 
production of UHE 7 —rays, and the neutrino source dis¬ 
tribution. Also, Ref. j 2 t| did not take into account prop¬ 
agation of UHE photons, instead assuming that all pho¬ 
tons are down-scattered into the GeV region. In addi¬ 
tion, simply due to the much larger statistics at lower 
energies, the quality of the fits performed in Ref. |2(| 
is dominated by the low-energy background component. 
Finally, Ref. j26| assumed that the sources do not emit 
any 7 —rays, although the 7 —ray energy fluence produced 
by pion production of accelerated nuclei should be com¬ 
parable to the pro duced neutrino fluence, as will be dis¬ 


cussed in Sec. Ill 


In the present paper we show that for all neutrino 
masses in the range 0.07eV < m v < leV one can find 
parameters that fit the UHECR observations with com¬ 
parable quality. We therefore conclude that, at least at 


the current state of knowledge, it is impossible to extract 
evidence for the relic neutrino background or even best fit 
values for absolute neutrino masses from UHECR data. 

We do not consider in the present paper neutrino in¬ 
teraction channels with multiple IV ± and/or Z° produc¬ 
tion. These channels could be important in case of neu¬ 
trino masses m > 3 eV jnj, which however are strongly 
disfavored by considerations on large scale structure for¬ 
mation | |28|| . 

By detailed numerical flux calculations we show that a 
non-uniform source distribution allows the Z-burst mech¬ 
anism to explain the UHECR flux without substantial 
relic neutrino over-densities. However, this only works if 
the sources exclusively emit neutrinos. Because isospin 
symmetry requires the energy fluence of neutrinos and 
7 —rays produced by hadronic charged primary inter¬ 
actions in the source are comparable, this will require 
the photons to be down-scattered below the GeV range 
within the source. 


II. NEUTRINO SOURCE 


We assume in this section that a pure neutrino source 
model can somehow be constructed and start with this 
case. In the next section we relax this condition and 
include other primary particles into consideration. 

Our simulations are based on two independent codes 
that have extensively been compared down to the level of 
individual interactions. Both of them are implicit trans¬ 
port codes that evolve the spectra of nucleons, 7 —rays, 
electrons, electron-, muon-, and tau-neutrinos, and their 
antiparticles along straight lines. Arbitrary injection 
spectra and redshift distributions can be specified for 
the sources and all relevant strong, electromagnetic, and 
weak interactions have been implemented. For details 
see Refs. |30|]3l| 1. Specifically relevant for neutrino in¬ 
teractions in the current problem are both the s-channel 
production of Z bosons and the t-channel production of 
W bosons. The decay products of the Z boson were 
taken from simulations with the |3^] Monte Carlo event 
generator using the tuned parameter set of the OPAL 
Collaboration j33|. The main ambiguities in propaga¬ 
tion concern the unknown rms magnetic field strength 
B which can influence the predicted 7 —ray spectra via 
synchrotron cooling of the electrons in the EM cascade, 
and the strength of the URB which influences pair pro¬ 
duction by UHE 7 —rays (34|] . Photon interactions in the 
GeV to TeV range are dominated by infrared and optical 
universal photon backgrounds (IR/O), for which we took 
the results of Ref. Q . The resulting photon flux in GeV 
range in not sensitive to details of the IR/O backgrounds. 

Predictions for the nucleon fluxes agree within tens of 
percents whereas photon fluxes agree only within a fac¬ 
tor ~ 2 between the two codes. The latter mostly reflects 
the ambiguities in photon propagation mentioned above, 
but has no influence on the conclusions of this paper. 
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For the present investigation we parameterize the neu¬ 
trino injection spectra per comoving volume in the fol¬ 
lowing way: 

ME, z) = /(I + z) m E~^Q{E v max - E) 

Zmin V Z < Zmax , (1) 

where / is the normalization that has to be fitted to the 
data. The free parameters are the spectral index q„, the 
maximal neutrino energy ax , the minimal and maxi¬ 
mal redshifts z m i n , z max , and the redshift evolution in¬ 
dex m. We assume for simplicity that all six neutrino 
species (three flavors including antiparticles) are com¬ 
pletely mixed as suggested by experiments j3(| and thus 
have equal fluxes given by Eq. (fD). Finally we chose the 
Hubble parameter Hq = 70 km s -1 Mpc ^ 1 and a cosmo¬ 
logical constant Ha = 0.7, as favored today. 



E[eV] 

FIG. 1. Fluxes of neutrinos, 7 —rays, and nucleons pre¬ 
dicted by the Z-burst mechanism for m, = 0.5 eV, assum¬ 
ing sources exclusively emitting neutrinos with fluxes equal 
for all flavors. We show three cases of the source evolu¬ 
tion parameter, m = —3,0,3 by solid, dashed, and dotted 
lines, respectively. Values assumed for the other parame¬ 
ters are: B = 10 ^ 9 G, minimal URB strength, z m i n = 0, 
Zmm = 3, = 2 x 10 22 eV, q v = 1. For each case the 

neutrino flux amplitude / is obtained from minimizing y 2 for 
E m i n = 2.5 x 10 19 eV. Also shown are experimental upper 
limits on 7 —ray and neutrino fluxes (see text and Ref. 0 for 
more details). 

For a given set of values for all these parameters we find 
the neutrino flux amplitude / in Eq. (|lj) obeying all ex¬ 
perimental bounds on photon and neutrino fluxes and ex¬ 
plaining the UHECR flux at highest energies above some 
value E m ; n by the secondary UHE protons and photons 
by a maximum likelihood fit. The fit quality is charac¬ 
terized by a y 2 value. Note that there are many differ¬ 
ent kinds of extragalactic sources which can contribute 
to the observed UHECR flux with energies below the 
GZK cutoff Egzk — 4 x 10 19 eV. Thus, one should take 
E m i n 1$ -Egzk if one wants to explain all UHECR data 
above the cutoff by the Z-burst model. 


Fig-1 illustrates how unrealistically high local neutrino 
background over-density could be avoided by assuming 
sources that are more abundant at low redshifts. In this 
figure we show primary neutrino and secondary proton 
and photon fluxes for the case m „ = 0.5 eV. The fol¬ 
lowing values have been assumed for the parameters of 
Eq. (1): spectral index q v = 1, maximal neutrino en¬ 
ergy E^ ax = 2 x 10 22 eV, minimal and maximal redshifts 
Zmin = 0 and z m ax — 3. Three cases, corresponding to 
the redshift evolution index m = —3, 0, 3 are plotted as 
solid, dashed and dotted lines, respectively. A typical 
value, B = 10 ~ 9 G, is assumed for the EGMF as well 
as the minimal strength consistent with observations for 
the poorly known URB |34). The latter results in op¬ 
timistic predictions for the UHE 7 —ray flux. The neu¬ 
trino flux amplitude was fitted as described above for 
E m in = 2.5 x 10 19 eV. Thus we require that the Z-burst 
model contributes to 9 bins of non-vanishing flux in the 
AGASA data. For all three cases in Fig. [j] we obtained 
X 2 (9) — 4. 

As one can see from Fig. |l], the value of the source 
evolution parameter m mainly affects photons with GeV 
energies. The value m = 3 which was chosen in previous 
work [^5| is similar to evolution of active galaxies. The 
secondary photon flux for such a source distribution is 
in conflict with the diffuse GeV photon background ob¬ 
served by the EGRET experiment. The uniform source 
distribution, m = 0 , is already in agreement with the 
EGRET flux, while negative values of the flux, m = —3, 
lead to GeV photon fluxes well below it. The latter case 
corresponds to sources which are more abundant now 
than at high redshifts. For example, BL Lacertae ob¬ 
jects for which correlations with UHECRs were found in 
Ref. [| 8 |, are distributed in such a way. Note that the 
choice of unknown EGMF strength B and the URB flux 
affect only the UHE flux of photons and is not important 
in the EGRET region which is only sensitive to the total 
injected EM energy. 

Fig. I shows results for varying neutrino masses with 
the unknown parameters B 1 radio background strength, 
as well as /, q„, E^ ax , z min , z max , and m chosen such 
as to produce fits to the UHECR data of comparable 
quality. We present in Fig. ||two extreme cases of small 
m v = 0.1 eV and high m v = 1 eV neutrino masses. 
Because parameter space is huge, we fix some parame¬ 
ters to given values (z m i n = 0, z max = 3, and minimal 
URB strength) and vary only the evolution parameter 
m, the EGMF strength B, and the maximum energy 
Emax f° r every given neutrino mass. Again, we deter¬ 
mine the neutrino flux amplitude / from minimizing y 2 . 
For m v = 0.1 eV we get y 2 ( 6 ) = 1.6, while for nh, = 1 eV 
we have y 2 (9) = 2.6, i.e. the fit qualities are compara¬ 
ble. For all intermediate masses we also find similar fit 
qualities. From this we conclude that no preferred values 
for the neutrino masses m v can be extracted. Instead, 
for every neutrino mass exists a large parameter region 
in which the Z-burst model with pure neutrino sources 
may work. 
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FIG. 2. Dependence of flux predictions in the Z-burst 
model on neutrino mass. (a) Small neutrino mass 
m, = 0.1 eV, with m = —3, B = 5x 10” n G, ax = 10 23 eV. 
(b) Large neutrino mass m, = 1 eV, with m = 0, 
B = 10 ~ 12 G, = 10 22 eV. See also Fig. § for expla¬ 

nations. 

In Fig. j|a we show the case of small neutrino mass, 
m v = 0.1 eV, with m = —3, B = 5 x 10 _ 11 G, and 
£JJ,a X = 10 23 eV. For small neutrino masses the reso¬ 
nance energy is large and thus secondary photons and 
protons are produced at higher energies, apart from the 
photons produced by t-channel leptons. Due to electro¬ 
magnetic cascades most of the EM energy ends up in the 
GeV region and thus the EGRET flux gives the most 
stringent bound. In particular, for m v <0.1 eV, the pa¬ 
rameter space for the Z-burst model shrinks, and even 
a m = — 3 distribution of sources is not consistent with 
the data. However, there is no pronounced cutoff for 
photons in this case (due to the a < 1 power law, see 
details in Ref. ©)• This allows to explain some fraction 
of the UHECR events by photons. Finally, the required 
neutrino flux is higher for small neutrino masses which 


makes the GLUE experiment bound on the UHE neu¬ 
trino flux J23|] a crucial constraint for the Z-burst model. 

In Fig. 9b we present the case of large neutrino mass, 
m„ = 1 eV, with to = 0, B = 10 _ 12 G, and E^ ax = 
10 22 eV. In this case the required flux of neutrinos is 
somewhat smaller. The available parameter space for 
the Z-burst mechanism is large and the EGRET bound 
on the GeV photon flux can be met even for a uniform 
distribution of sources. 


III. PHOTON AND NEUTRINO SOURCE 


It is well known that sources capable of accelerat¬ 
ing UHECRs produce 7 —rays up to at least U^ ax ~ 
100 TeV j37j . Since in acceleration scenarios both 7 —rays 
and neutrinos are produced as secondaries, the power 
/ 7 ~ 1 radiated in 7 —rays relative to neutrinos has to 
be of order unity. For the 7 —ray injection spectrum this 
implies 

ME, z) = 6// 7 (1 + E 0(E^ ax - E) 

9\Qrf) 

^min ’max • (2) 

Here, q 1 is the 7 —ray spectral index, and g(q, E max ) = 
f5 max dE U 1-9 , where we have introduced some small 
low energy cut-off U m i n for convergence. 



n E [eV] 

FIG. 3. Same as Fig. Ml (m„ = 0.5 eV) for sources emitting 
equal power in neutrinos and 7 —rays, / 7 = 1 in Eq. (Q). 
The other parameters are chosen as m = 0, q-y = 2, 
Ffkax = 100 TeV. 

In this section we consider the more realistic case 
where the sources also emit 7 — rays with a power com¬ 
parable to the emitted neutrino power, / 7 = 1, up to 
EJax = 100 TeV in Eq. ([|). For all other parameters we 
chose values that minimize the tension with the obser¬ 
vational upper limits on the UHE neutrino flux and the 
diffuse GeV 7 — ray flux. Fig. [| shows that in this case the 
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Z-burst scenario cannot be made consistent with obser¬ 
vations. A possible solution to this problem is to down- 
scatter most of the EM energy into the sub-MeV range 
within the source. Only in this case can the EGRET 
bound be satisfied. This would requires a very strong 
photon field up to >keV within the source. 

Even the scenario in Fig. |] is still optimistic because it 
assumes that the source is completely opaque to the pri¬ 
mary nucleons. While this may be achieved easier than 
containment of 7 —rays, for example, by magnetic fields, 
it is clear from Fig. |3j that even if only a small fraction 
of the primary nucleons leave the source, the nucleon 
flux between ~ 10 18 eV and ~ 10 19 eV would be much 
higher than observed, in agreement with the conclusions 
of Ref. Q . This problem could be avoided if the pro¬ 
tons are deflected strongly enough so that they could not 
reach the Earth. However, this possibility also appears 
unrealistic as has been discussed in Ref. |35J]. This prob¬ 
lem can only be solved if the protons are trapped in the 
source. 

Finally we note that the Z-burst mechanism also poses 
extreme requirements on the acceleration mechanism it¬ 
self since the primary protons have to be accelerated 
to energies -ET ax ~ 10E^ lax > 4 x 10 22 eV (eV/m u ). 
In contrast, known mechanisms are usually limited to 
^ ax < 10 22 eV@. 

Thus, the Z-burst model imposes the following require¬ 
ments onto the sources: They should emit energy only in 
neutrinos and in sub-MeV 7 —rays, and should also trap 
most of the primary protons. 


burst mechanism would exhibit a GZK-cutoff for nucle¬ 
ons and would be dominated by 7 —rays at higher en¬ 
ergies. Furthermore, the required UHE neutrino fluxes 
are close to existing upper limits and should be eas¬ 
ily detectable by future experiments such as Auger [ (fl] , 
Euso fu|, RICE JiH , or by other radio detection tech¬ 
niques [Q . 

The space of parameters characterizing neutrino 
sources and their evolution is highly degenerate when 
fluxes are fit to the observed UHECR fluxes. Since 
evidence of relic neutrinos and extraction of absolute 
neutrino masses requires conservative assumptions about 
these unknown parameters, we conclude that the current 
state of knowledge does not allow to extract any meaning¬ 
ful information on neutrino masses from UHECR data. 
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IV. CONCLUSIONS 


The Z-burst mechanism where the highest energy cos¬ 
mic rays are produced by neutrino beams interacting 
with the relic neutrino background only works with 
sources exclusively emitting neutrinos in the ultra-high 
energy regime. In order to avoid conflict with the known 
diffuse backgrounds of 7 —rays, these sources should emit 
photons only in the sub-MeV region. In addition, they 
should trap primary protons in order to avoid an exces¬ 
sive nucleon flux from the source, and should be able to 
accelerate these protons up to £T ax > 10 23 eV (eV/m„). 
None of the astrophysical acceleration models existing in 
the literature seems to meet this requirement. 

Under the assumption that such an extreme source 
class nevertheless exists we have shown that the Z-burst 
mechanism can work without unrealistically high local 
relic neutrino over-densities if the neutrino sources are 
typically more abundant at present than in the past. 
Especially neutrino masses m v < 0.5 eV require a non- 
uniform source distribution oc (1 + z) m with negative 
evolution factor, to < 0, as is the case with BL Lacertae 
objects. 

The contribution to the UHECR flux from such a spec¬ 
ulative extragalactic neutrino source class due to the Z- 
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